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A method is presented which allows the exact construction of conserved (i.e. divergence-free) 
current vectors from appropriate sets of multipole moments. Physically, such objects may be taken 
to represent the flux of particles or electric charge inside some classical extended body. Several 
applications are discussed. In particular, it is shown how to easily write down the class of all smooth 
and spatially-bounded currents with a given total charge. This implicitly provides restrictions on 
the moments arising from the smoothness of physically reasonable vector fields. We also show that 
requiring all of the moments to be constant in an appropriate sense is often impossible; likely limiting 
the applicability of the Ehlers-Rudolph-Dixon notion of quasirigid motion. A simple condition is 
also derived that allows currents to exist in two different spacetimes with identical sets of multipole 
moments (in a natural sense). 
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I. INTRODUCTION 

There are a number of instances of "currents" appear- 
ing throughout physics as fields which somehow express 
a conservation law. Here, we will be studying vector 
fields 3(x) in 3 + 1 dimensional spacetimes. In classical 
physics, a common example of such an object would be 
the electromagnetic current. Another would be the nu- 
merical flux density of some species in a continuum (as 
long as the total number of such particles isn't expected 
to change via chemical or other reactions). 

Keeping these examples in mind, we take the sup- 
port of such currents to be a spatially-bounded region W 
which surrounds a timelike worldlinc Z. Intuitively, this 
could be thought of as the worldtube of some extended 
body. The "central" worldline would then be most nat- 
urally identified with the body's center-of-mass. 

We may also use any spacelike hypersurface S which 
cuts W to define a charge in the usual way: 

q(S) := f J dS . (1.1) 
Js 

The current would be "conserved" if q{S) were indepen- 
dent of S. This clearly occurs when 

V-J = 0. (1.2) 

Our purpose in this paper is to study vector fields 
which satisfy this equation and have the given support 
properties. In a sense, finding such fields is actually quite 
simple. This can be seen by constructing a 3-form from J 
in the natural way, and then observing that (|1.2|) is equiv- 
alent to requiring that this form be closed. Solutions fol- 
low immediately. For example, [l[ discusses a common 
procedure in continuum mechanics whereby the number 
density flux is essentially the pullback to spacetime of a 
fixed volume form defined in a three-dimensional "body 
space." This method is quite natural when discussing the 
worldlines of a body's constituent particles, but may be 
rather awkward in other cases. 



For example, it is often useful to decompose a current 
into a set of multipole moments. In some cases, mo- 
ments may be experimentally determined be measuring 
the asymptotic behavior of the fields sourced a body. If J 
happened to be an electric current, the moments would 
also appear as coefficients in an expansion for the net 
force and torque due to externally applied electromag- 
netic fields S| . Even when including the effects of 
self-fields, some aspects of the multipole decomposition 
remain useful in a wide range of scenarios 0, [E[ • 

These types of methods also provide an elegant way of 
separating out different scales when discussing the struc- 
ture of some extended body. Since the moments (that we 
will choose) are in fact tensor fields defined along Z, it 
becomes relatively simple to discuss this structure from 
the point of view of a single comoving observer. This 
makes a number of problems more intuitive. 

Unfortunately, multipole moments are related to J (or 
the deformation gradient which may be used to construct 
it) in a rather complicated way. In this paper, we shall 
take the point of view that the moments themselves are 
fundamental. Other more commonly-discussed quanti- 
ties are to be derived from them. It is shown how to 
explicitly construct J in this way in Sec. |TT1 which leads 
to the introduction of two "moment potentials" defined 
on the tangent bundle restricted to Z. A moment may 
be derived from these quantities essentially by taking a 
(three-dimensional) Fourier transform, and then differen- 
tiating a sufficient number of times. Importantly, these 
potentials are not required to obey any evolution equa- 
tions. Their "time-dependence" is essentially arbitrary. 

Not all of them will describe physically reasonable cur- 
rents, however. Most choices introduce singularities. Sec. 
IIIII discusses how to avoid such problems. This implicitly 
shows which restrictions must be placed on the moments 
for them to represent plausible current vectors. 

As a byproduct of this construction, we obtain a sim- 
ple method of constructing the class of essentially all 
smooth currents with a given total charge. The situation 
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is nearly identical to the one in flat spacetime, where the 
charge density and three-current observed by someone on 
Z are very simply related to the moment potentials. This 
could be useful in studying corrections to the Dewitt- 
Brehme expression |6j for the electromagnetic self-force 
on a charged test particle in curved spacetime (which is 
intended to describe all "sufficiently small" charges with 
a given q). The flat spacetime analog of this idea was 
carried out in 

Lastly, we investigate some consequences of imposing 
simple a priori conditions on the moment set itself. Sec. 
IIVI supposes that the moments are all constant with re- 
spect to some frame attached to the central worldlinc 
Z. This closely mirrors the idea of quasirigid motion ad- 
vanced in @, 0, Is! ■ Unfortunately, it is shown that the 
moments can only remain constant under certain condi- 
tions involving the Riemann tensor and the acceleration 
of Z (roughly speaking). 

Sec. IVlthen looks at the possibility of two currents - 
possibly in different spacetimes - having identical sets of 
multipolc moments as viewed in their respective frames. 
As in the case of setting all of the moments constant, 
this is not always possible. A restriction on the Riemann 
tensors on both central worldlines is derived which allows 
these types of "comparable" currents to exist. 

The system of multipole moments adopted here is the 
one developed by W. G. Dixon |, 0, 0, i, ■ The rele- 
vant portions of the formalism are reviewed in appendix 
|B1 It has a number of desireable properties without which 
many of the results in this paper would likely be impos- 
sible to reproduce. It also has a very elegant place in the 
(exact) theory of relativistic mechanics. And perhaps 
most importantly, the formalism also works for symmet- 
ric rank-2 tensors T satisfying V • T = (or more gen- 
erally V • T = F • J with F a closed 2-form; e.g. an elec- 
tromagnetic field) . This paper may therefore be taken as 
a model for attempting to construct all possible stress- 
energy tensors with a given center-of-mass and net linear 
and angular momentum vectors. 

The results here also make some use of bitensors. Ap- 
pendix [A] reviews the subject, and states some results 
used in text. 

For the remainder of this paper, indices shall always be 
included on tensor or bitensor quantities. Those which 
refer to points on Z will be denoted by a, b, .... Indices 
associated with all other points in W will be written as 
a', b', .... Tetrad labels running from to 3 are to be 
expressed in capital Roman letters from the beginning of 
the alphabet, while triad labels running from 1 to 3 will 
be denoted by I, J, . ... We choose the signature of the 
metric to be — h ++. 



II. CONSTRUCTING CURRENTS 

At this point, we assume that the reader has at least 
skimmed the appendices to gain familiarity with the no- 
tation and basic elements of Dixon's formalism. These 



outline how conserved vector fields in spacetime may be 
represented by the current skeleton J a ; a vector- valued 
distribution defined on the tangent bundle (restricted to 
a neighborhood of Z). Dixon has actually proven that 
all well-behaved currents may be uniquely represented 
in this way [H . 

The reverse is not true, however. One could pick mo- 
ments satisfying all of the given constraints, and yet the 
right-hand side of (|B9[) may not even exist. Then there 
wouldn't be any current vector at all associated with that 
set of moments. Even when the integrals do exist, the re- 
sulting J a may be singular. Such a result would clearly 
be unphysical. 

The main purpose of this paper is to remove this prob- 
lem. We want to be able to specify a current skeleton 
in some simple way which automatically guarantees that 
the associated vector field is both physically reasonable 
and has a given charge q. We do this by constructing 
J a explicitly in terms of J a , and then checking that the 
result is reasonable. 

Our first step will be to construct current vectors from 
all of the skeletons satisfying the constraints described 
in appendix [Bj The natural starting point for this is 
of course (|B9[) . It is convenient to convert the integral 
over X on the right-hand of this equation into one over 
x = exp z X. The Jacobian of this transformation is given 
by (|A12[) . from which it follows that 

(J,0) = J ds(Aj a H a ' a ,<f> c ,) . (2.1) 

As in (|B9[) , <f> a i is an arbitrary test function. This equa- 
tion therefore determines the current completely in terms 
of its skeleton (ignoring irrelevant sets of measure zero). 
Assuming that the integrals on the right-hand side are 
absolutely convergent, they may be freely commuted. 
Taking advantage of this, one finds that 

J a ' - j dsAJ a H a ' ,^ A =0 , (2.2) 

which allows the identification 

J a \x) = J dsA(x,z)J a (X 7 z)H a ' a (x 1 z) . (2.3) 

Here, X a = —a a (x, z) and z = z(s). 

In [2| , it was found that the constraint equations (|B3[) - 
(|B5|) imply that J a - defined in the text immediately 
preceding (|B16[) - has the form 

J a (X,z) = A a {X,z)5(n ■ X) - h a b (z)X b 

xB(X,z)6'(n-X) . (2.4) 

h a t, := S^ + n a rib is the standard projection operator, and 
it is given that n b V \bA a — n h V*b£> = 0. The only other 
restrictions on the moments may be written as 

n a A a = V w (h a b X b B) , (2.5) 
V« a A a = . (2.6) 
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Since A a does not depend onn ■ X, this second equation 
is really only a 3-divergence. Also note that by definition, 



J a = Q a S{X) + j a 



(2.7) 



Using these expressions, the s-integral in (|2.3j) can be 
performed explicitly to yield 



J a ' = / dsS(X)Q a ' +N- 1 



AA a H a 



ABh a b X b H a a 



(2.8) 



where the right-hand side is understood to be evaluated 
at the (unique) z which sets n(z) ■ X(x, z) — 0. Also, 



N 



-n a a a b v b + h a X a 



(2.9) 



Following the notation of (|B1[) . h a := 5n a /ds. In flat 
spacetime, this expression reduces to the lapse of the foli- 
ation {£}. Although it doesn't retain this interpretation 
more generally, we shall still refer to it as a lapse. 

Now for x £ Z, note that the integral in (|2.8p doesn't 
contribute to the current. Temporarily restricting our- 
selves to such cases, (|A5|) and (1A13[) can be used to show 
that 

J a \x) = AN^H 11 ' a \A a + X a v b V b (B/N) 



(B/N) (n a h b X b - h a b a b c v c 



-2H b ' h v c xK b l J> , 



(2.10) 



It is natural to split this up into components adapted to 
an observer moving along Z. Define the charge density 
to be 



p(x) := n a (z(r))cr a a > (x, z(r))J a (x) . 
Similarly, let the three-current be given by 
f(x) := -h a b (z(T))a b b ,(x,z(T))J b '(x) 



(2.11) 



(2.12) 



Substituting into ([2TTD]) . and using (j23|) and (j!H)D . the 

charge density reduces to 



-AV* n [h\X b (B/N)] 



(2.13) 



With the exception of the overall factor of A, this is ex- 
actly the same form that occurs in flat spacetime. Such a 
simple correspondence will not hold for the three-current, 
however. 

Before computing it, note that B does not appear 
by itself in (|2.13p . One instead finds the ratio B/N. 
This grouping will continue to be useful when expand- 
ing the three-current, although it is a little inconvenient 
as stated. We would like this ratio to have the same argu- 
ments as B itself. Specifically, it should be independent 
of n ■ X. But n a \7* a N(X, z) ^ even when n ■ X = 0. 
To fix this, define a "flattened" lapse N as 



N(X,z) :=N(h{z) -X,z) , 



(2.14) 



and a "normalized" moment potential 

C{X,z) := B(X,z)/N(X,z) 



(2.15) 



Clearly, n a V* a C(A, z) — 0, as desired. The projected 
divergence in (12.13[) allows C to be substituted for B/N 
there, so 



-AV* a (h a b X b C) 



(2.16) 



This may be further simplified by introducing an or- 
thonormal tetrad {n a (z) , ef (z)} defined on Z. Enforcing 
the orthonormality conditions requires that the spatial 
triad evolve according to 



e"j := v b \7 b e a j = n a h b e b j + fi/e} , 



(2.17) 



where Qu(z) = Q,nn(z) is arbitrary, and represents the 
angular velocity of the triad. It needn't have any par- 
ticular physical significance, however, so we shall leave it 
unspecified for now. 

If the tetrad is assumed known, we can trivially intro- 
duce a set of (normalized) Riemann normal coordinates 
defined by 



R A (x,z) := e£ (z)X a (x,z) 



(2.18) 



for any fixed z. Then an arbitrary bitensor f(x,z) could 
be rewritten as either (abusing the notation a bit) f(X, z) 
or f(R,z). 

This is particularly useful for objects such as C, which 
are independent of R° = —n ■ X. So we may write 
C = C(R, z), where R represents the spatial coordinates 
{R 1 }- Defining the differential operator di in the obvious 
way, p,16p now simplifies to p = -A<9j (i? J C). 

Recall that this expression is only correct away from 
Z. Correcting it requires evaluating the integral in (|2.8[) . 
This is easily done by noting that 6(X) = 6{R°)5 3 {R). 
Then 



dsS(X)Q a =5 3 {R,)Q a 



(2.19) 



Using (|B6|) , the charge density everywhere is finally given 
by 



p = A [q5 3 (R) - d z (R*C)] 



(2.20) 



It is now useful to define a function ^(R, z) which satis- 
fies 



di = 5 3 (R) 



(2.21) 



This equation will be solved in the following section. But 
for now, tp will be left unspecified. 
We instead define 



tp(R, z) := #(R, z) - C(R, z) 



(2.22) 



Anywhere in a neighborhood of W , the charge density is 
now given by 



p = Ad! 



(2.23) 
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B (and therefore C) was originally defined to vanish 
outside of W. But applying Gauss' theorem to (|2.2ip 
shows that ip cannot share this property, tp must there- 
fore be equal to qip everywhere outside of W. This clearly 
implies that p — in such regions, as required. Beyond 
this restriction, ip may be specified arbitrarily as long as 
the charge density which it determines remains nonsin- 
gular. 

Implicit in these statements is the fact that any func- 
tion depending on (EL, z) may be uniquely identified with 
another function depending only on x. In a sense, we are 
simply choosing to work in a special system of coordi- 
nates. These happen to be closely related to Fermi nor- 
mal coordinates, and exactly coincide with them when 
n a = v a and Q IJ = 0. 

In any case, we may now apply a similar analysis to the 
three-current. Using (f2T8|) . (|2T0|) . and (|2~T2|) . the tetrad 
components of j a reduce to 



= AN- 



A 1 + g?/<5 3 (R) + R T C + dj (i? J C) 



(n J K R K - eta b c v c ) -C^jR 3 
2e c jH c ' c v d a b dc ,e [ *R^ 



(2.24) 



This uses the notation j 1 := e^j a , A 1 := e I a A a , etc. But 
it can be rewritten more suggestively using (|2.20[) : 



AJV- 1 !^ + R J C - (p/A) {rtjR 1 



e{v b a a b ) 



2d, 



R [I (n J] K R K - e$v b (T a b y~\ } . (2.25) 



This form is useful because the second line is manifestly 
(three-) divergence- free. From (|2.6p , A 1 also has this 
property, but is otherwise arbitrary (ignoring its support 
properties for the moment). We may therefore absorb 
these terms by defining 



K 



J 1 := A 1 + qR 1 ^ + 2dj R lI (n J] K j 

-e J Jv b a a b )c] . (2.26) 



The three-current then reduces to 



J 1 = N-^AiJ 1 -R'p) +p(e I a v b a\ 



rtjR J 



(2.27) 



In summary, a conserved current vector may be con- 
structed by choosing "moment potentials" <p(R, z) and 
J r/ (R, z). ip must reduce to qip outside of W, but is oth- 
erwise arbitrary if we don't worry about smoothness of 
J a . Similarly, J 1 must reduce to qR Ir ip outside of the 
worldtube. It is also required to satisfy diJ 1 = every- 
where. Given functions with these properties, the current 
is straightforward to compute from (|2.23p , (|2.27|) , and 



J a ' =H a ' .(n°p + e?/) 



(2.28) 



This last equation is a direct consequence of (|2.1ip and 

§21. 



III. CHOOSING if AND J 1 

We saw in the previous section that nearly any current 
may be constructed from potentials ip and J 1 chosen to 
satisfy very simple algebraic matching conditions. But 
many such choices generate singular current vectors. In 
this section, we show how to exclude this class; allowing 
one to easily construct physically reasonable divergence- 
free vector fields. 

Since ip must reduce to qip outside of W, our first step is 
to explicitly solve (|2.2ip . Adopting spherical coordinates 
(r, 9, <f>) defined in terms of R in the standard way, the 
solutions of this equation are easily seen to be of the form 



1 



47rr 3 



1 



■q-iAfafrz)) 



(3.1) 



where A(9 1 </>; z) is any function with a vanishing 
monopole moment (in the elementary sense). The fac- 
tor q~ x is for later convenience. 

The arbitrariness of A simply means that ip is not 
uniquely determined by (|2.2ip . which is to be expected. 
This ambiguity translates into a freedom in ip as well. 
The same conclusion may also be reached directly from 
(|2.23p , where a given p only determines ip up to terms of 
the form A/r 3 . But there cannot be a 1/r 3 divergence 
in tp (without any angular dependence), as this would 
imply that the charge density contained an unphysical 5- 
function at the origin. It follows that there always exists 
a ip which completely removes any o(r~ 3 ) singularities 
in ip. From now on, we always choose potentials in this 
class. This choice actually removes the only polynomial 
singularity that may exist in tp. Others would lead to 
unphysical charge densities. 

Most (^-function singularities are excluded from ip on 
similar grounds. The only exceptions are terms of the 
form a(z)6 3 (R). But these do not affect the current or 
the moments in any way, so they have no physical rele- 
vance. We can therefore set a(z) = without any loss of 
generality. 

Now choose a ip that is continuous in both R and z. 
As long as the boundary of the worldtube can be writ- 
ten in the form r = 71(9, <p; z), where 1Z is single- valued, 
specifying ip in some neighborhood of W determines ip 
everywhere: 



(3.2) 



The problem with this approach is that arbitrary choices 
of ip will often lead to charges that vary over time. 

It is much simpler if we look for a ip that fits a given 
ip rather than a ip which fits a given ip. This has the 
advantage of allowing one to specify the charge a priori. 
It is then trivial to ensure that it remains constant. 

Before showing how to do this, we first require that the 
charge density at r = does not depend on the angular 
variables (i.e. that it is single- valued) . Define a central 
density po(z) such that 



p (z) := p(0,9,(p;z) . 



(3.3) 
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Now in order to find <p, choose q, A, p , and TZ(6, (f>; z). 
For r < TZ, <p may then be taken to have the form 



'•P 



1 r 



(3.4) 



(p(r, 9, </>; z) is any bounded and continuous function in 
W. As an example, a charge density of poA is generated 
by choosing (p = and 



47T , 



(3.5) 



One sees that A vanishes in this case when the boundary 
is a sphere. 

With very little extra effort, we could also take into 
account the charge density on the body's surface, or any 
number of other boundary conditions. This isn't particu- 
larly necessary, however, and the extensions are straight- 
forward. 

The remainder of the current vector is now found 
by specifying an arbitrary divergence-free "three- vector" 
J r/ (R, z) which reduces to qR 1 ^ outside of the body. 
If we choose a nonrotating triad, this vector field has 
a fairly straightforward physical interpretation in many 
cases of physical interest. Suppose that n a and v a are 
determined by center-of-mass conditions, and that the 
forces and torques on the body are not too large. Then 
n a ~ v a , and the components of n 1 would be small rela- 
tive to the (reciprocal of the) worldtube's proper diame- 
ter. If the spacetimc curvature inside W is also small in 
an appropriate sense, (|2.27|) will reduce to j 1 ~ J 1 —R 1 ^. 
So physically, J 1 may be interpreted as essentially that 
portion of the 3-current which arises due to the bulk 
motion "internal" charges. In this approximation, J 1 
should therefore be nonsingular throughout the body. A 
slightly more detailed look at (|2.27| and (|3.4| shows that 
this observation must remain true even in the exact the- 
ory. 

We are therefore left with the problem of writing down 
a nonsingular divergence-free three-vector which has a 
specified form outside of a given boundary. Finding vec- 
tor fields with these properties may be nontrivial when 
A / 0, so we now illustrate a method that may simplify 
the process. First write down J 1 so that it automatically 
satisfies the appropriate boundary condition: 



J 1 = Q(R - r)jfa + 0(r - TZ)qR 1 ^ 



0(x) 



(3.6) 
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Here, O(-) is the Heaviside step function; 
x > 0, but otherwise vanishes. 

This expression also introduces the new vector field 
(r, 9, </>; z) which must be both nonsingular and 
divergence-free. While we could require this to exactly 
match qR Ir ip at r = TZ, doing so would eliminate the pos- 
sibility of currents flowing on the surface of the body. 
Moreover, the divergence of (|3.6p will vanish everywhere 



if only the normal component of J 1 is continuous across 
the boundary. We therefore suppose only that 



A 



■R 



AttR 2 



where 



SI := d l (r - U) 



(3.7) 



(3.8) 



The boundary condition (J3T7J) is difficult to apply di- 
rectly if we (say) write t7/_\ m terms of a vector potential. 
To fix this, we instead construct it from two scalar po- 
tentials Q and P in such a way that its divergence will 
always vanish. To simplify the notation, it is first useful 
to define the operator 



e IJK R J d K 



(3.9) 



or "L :=Rx V" in a more traditional notation. This is 
just (i/H multiplied by) the standard angular momentum 
operator used in quantum mechanics. 
We may now write 



(-) 



JJK 



djh 



Ki- 



ll P 



(3.10) 



This expansion is essentially an application of the angular 
momentum Helmholtz theorem [ll| . The scalars may 
also be interpreted as Debye potentials [H, EH > or t ne 
generators of the toroidal and poloidal components of 

The utility of this expansion can be seen by considering 
scalar potentials which may be split into angular and 
(normalized) radial components. For example, let 



Q(r,6, (f>; z) 



\Qe{r/n;z)Qt{e,<f>;z) , (3.11) 



and break up P in the same way. We choose the first term 
to generate the right-hand side of (|3.7p . For £ > 1, the 
angular functions Qi could then be taken as appropriate 
basis functions. The corresponding Qg will be chosen to 
ensure that the £th term in the sum doesn't add anything 
to the normal component of Jf_\ at the boundary. 
Looking at the £ = terms, p.7[) becomes 



A_ 

4tt 



= L • 



QoKLQ - Q' Q LTZ 



-P P TZ 2 VTZ , 



(3.12) 



where Q' indicates a derivative of Qo with respect to 
r/TZ. Also note that these radial functions are under- 
stood to already be evaluated at r/TZ = 1 (i.e. L docs 
not act on them here). 

The solutions to (|3.12|) are clearly not unique (even 
taking into account the gauge freedom). Still, we only 
need to generate a solution to this equation. Supposing 
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that P =P = 0, Q (l;z) = 1, and Q' (l; z) = -I, we 
are left with the relatively simple result 

^ = L 2 (kQo) . (3.13) 

This may be solved by expanding IZQo in spherical har- 
monics, using the addition theorem to express the series 
in terms of Legendre polynomials, and then summing the 
result. One eventually finds that (up to a possible addi- 
tive function proportional to 1 /1Z) 

Qo{9^-z) = -—^ J dft'i(#'>';z) 

x ln(l - cos a) . (3.14) 

df2' is that standard volume clement on a sphere, and 
a is the angle between the point of integration and the 
point where Qo is being evaluated. 

The radial profile Qo may now be chosen. This must 
satisfy the aforementioned conditions Qo(l; z) = 1 and 
Qo(l; z) — —1, but we must also require 

\im Q (r/K; z)/r = . (3.15) 

r^O 

This is necessary to ensure that the three-current is non- 
singular and single-valued at the origin. A particularly 
simple choice is then 

Qo(r/K;z) = (£) 2 ( 4 -3^) . (3.16) 

This completes our solution to (|3.12[) . 

But as was already mentioned, the full solution to (|3.7[) 
is usually more complicated. If each term in (I3.11|) may 
be considered individually, then for all £ > 1, 

QiKUQi = Q' e QiUK + P e PiK 2 diK (3.17) 

at r = 1Z. 

If dilZ = 0, this implies that either Qe(l;z) vanishes, 
or that Qe is independent of 9 and <j). But throughout the 
body, only the angular derivatives of Qe have any phys- 
ical relevance, so this case is equivalent to Qe = 0. But 
this is itself indistinguishable from saying that Q vanishes 
everywhere. So we may always say that Qe vanishes at 
r = 1Z when the boundary is locally spherical. 

In this case, J(_^ reduces to a sum of terms of the form 

Qe (jj^^Qi ~ d'Qe^j + ^Q'e^Qe + Pe^ 1 Pi ■ (3.18) 

On the boundary, the first group of terms clearly vanishes 
foralH>l. 

But we must also make sure that f|3 . 18[) is everywhere 
finite. We first require none of the functions Pe contain 
any singularities. Choosing Q is slightly more compli- 
cated due to the o(r~ 1 ) divergence in (J3TTSJ). But for all 



I ^ 0, we may suppose that there exists a nonsingular 
function Q Y e such that 

Qiir/Tl; z) = ^ (l - ^) Qr,i{r/K; z) . (3.19) 

This guarantees that (|3. 1 8[) remains finite. As long as the 
sum over t does not diverge, these choices are sufficient 
to describe a wide variety of functions J 1 (for angles) 
where d{R. = 0. 

If the boundary radius is changing, the situation is 
more complicated, and our discussion will not be quite 
so complete. The full equation (|3.17p must be solved in 
general at r — 7Z, although it isn't obvious how to do 
this. To make some progress, suppose that Pe vanishes 
on the boundary. If Qe(l;z) also vanishes, Qe may be 
anything at all (so long as it remains finite). But more 
generically, it must have the form 

Qe oc^[^ (1;z)/c5 ' (1;2) ] . (3.20) 

The general case where Pe ^ at r = 1Z does not 
seem simple to discuss. It is likely that an ansatz better 
adapted to nonspherical boundaries than (|3.10[) would be 
useful for this purpose. The class of currents which we 
are excluding by neglecting this case does not really seem 
significant, however. 

To summarize, J 1 is required to compute the three- 
current. It has the form (|3.6p . where J}_\ may be ex- 
panded using (|3.10[) . The scalar potentials in this expres- 
sion have the form (|3 . . By definition, Pq = 0, and Qo 
and Qo are given by (|3.14|) and (|3. 16[) respectively. The 
higher-^? components of the potentials may be chosen by 
first writing down any (smooth) sets of functions Pe, Pe, 
and Qe. Then each Qe should have the form (|3.19[) . Note 
that while these choices do not allow for every possibility, 
they are still quite general. 

In all of this, the z-dependence of ip and J 1 is essen- 
tially arbitrary. Since z(s) is a one-to-one function, these 
potentials may easily be rewritten as functions of (R, s) 
rather than (R, z). When doing so, it seems sufficient 
to require that they be C 1 in s. All of these conditions 
taken together allow us to easily generate a very wide 
class of physically reasonable current vectors with a given 
q. Note that the continuity restrictions are sufficient but 
not necessary. Still, the stated classes seem large enough 
for almost any practical purposes. 

Now that acceptable classes of ip and J 1 have been 
identified, expressions in Sec. |TT] and appendix [B] may 
be used to find the which sets of moments correspond 
to plausible current vectors. Since translating between 
these two representations involves a Fourier transform, it 
is likely to be difficult to do this explicitly. We shall not 
attempt it, but merely state that most of the restrictions 
on the moments would involve their asymptotic behavior 
in the "large-n" limit. 
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IV. DYNAMICAL RIGIDITY 

In Newtonian mechanics, one often makes use of the 
concept of rigid bodies. Although such objects do not 
(and cannot) exist, the idea is still useful. In many cases, 
the bulk dynamics of a large class of real materials are 
known to be adequately approximated by the behavior of 
a suitable rigid body. Just as importantly, rigidity con- 
ditions allow equations of motion to be solved uniquely 
in relatively simple ways. Although it is not at all clear 
that a similarly useful concept should exist in highly rel- 
ativistic systems, it is still interesting to search for one 
that is roughly analogous. 

Various definitions of rigidity in curved spacetimes 
have been introduced over the years. But perhaps the 
simplest is a generalization of the concept advanced by 
Born for objects in Minkowski space [TH, [HI]. This 
essentially requires that the strain tensor throughout a 
material be Lie-dragged along the worldlines of its con- 
stituent particles. Unfortunately, such a restriction fails 
to be self-consistent in many systems of practical interest 

m- 

Rather than defining rigidity through local kinematics, 
it has been conjectured that fixing certain "global" (or 
quasilocal) quantities may be more fruitful. In partic- 
ular, one may wish to call a body rigid if its multipole 
moments evaluated in some particularly natural frame do 
not change. 

Supposing again that we are discussing the mechanics 
of some extended body, the reference line Z and the dy- 
namical velocity n a of this frame may be chosen using the 
center-of-mass conditions given in [§|, Q (and rigorously 
justified in [Hj]). In doing so, the temporal component of 
the tetrad would be proportional to the body's bulk (lin- 
ear) momentum vector. The triad vectors would not be 
chosen uniquely in this scheme, although their evolution 
may be fixed by defining an "average angular velocity" 
as the vector which must multiply the "inertia tensor" to 
yield the bulk angular momentum vector 0, Q . This ef- 
fectively defines fl 1 j in terms of the body's stress-energy 
tensor. 

These ideas were introduced in |9j under the name of 
quasirigidity, and supposed that the multipole moments 
under discussion were those of the object's stress-energy 
tensor. Very similar ideas have also been referred to as 
dynamical rigidity in @, H[, and we shall use the two 
terms interchangeably. These latter papers also intro- 
duced a scalar defined from the body's stress-energy ten- 
sor (and the metric) which may be interpreted as its total 
internal energy. It was found that this is in fact a con- 
stant of motion for bodies which are dynamically rigid. 
The concept therefore has an elegant physical interpre- 
tation. 

Changing definitions slightly, the results of the previ- 
ous sections may be used to explicitly construct quasi- 
rigid bodies. We do so by requiring that the multipole 
moments of the number density and electric current vec- 
tors associated with a given object be fixed in the afore- 



mentioned frame (rather than the moments of its stress- 
energy tensor). This is actually far less restrictive than 
the definition in d, H, as it ignores the stress and en- 
ergy density distributions inside the body. In doing so, it 
also appears closer in spirit to the concept of Newtonian 
rigidity. A downside is that the total internal energy of 
the body will no longer be conserved in general. It is also 
somewhat less "dynamical" than the original definition. 

There is one potential ambiguity in this: it might be 
useful to allow the monopole moment to vary even if the 
higher moments do not. If Q A := e A Q a were held con- 
stant, (|B6|) implies that qd(v I )/ds — 0, or 

q (n 1 - v 1 - Q 1 jv J ) = . (4.1) 

This is of course satisfied trivially if either q = or n a = 
v a and $Vj = 0. But in general, it provides a rather 
strict restriction on the class of allowed frames. 

In a mechanics problem, one would usually derive v a 
from n a using the mass, force, torque, and linear and 
angular momenta of the body [1,0. In these cases, (14. ip 
would restrict which combinations of these objects could 
be allowed, which seems unphysical. In any case, one of 
the main reasons for adopting rigidity conditions in the 
first place is to obtain unique evolution equations for all 
relevant quantities. But changes in Q a would already 
be fixed in most cases of interest. There is no need to 
force a potentially contradictory condition on top of the 
already existing one. We therefore define a current to be 
quasirigid if all of its moments except for (possibly) the 
monopole are constant. 

It is clear from (|B2[) that this definition implies that 
J a {X, z)e A {z) — J A is a function purely of e A X a = R A . 
Noting (|2.4p . this means that A 1 and B must be inde- 
pendent of z (or equivalently, of s). 

We first consider the effect of holding B constant. Sup- 
posing that p is at least instantaneously smooth, (|2.15[) 
and (|2.20p can be used to show that 

p~A(fl/Jv)^(iJ z 0/ln#) , (4.2) 

where terms that are obviously nonsingular have been 
suppressed. 

In general, B diverges at the origin. The charge density 
may therefore develop a singularity of its own if we are 
not careful. One necessary (but not sufficient) condition 
that must be satisfied to avoid this is 




Writing this in a more useful form requires an expan- 
sion of N in powers of R or X around r = 0. But that 
in turn requires a similar expansion for er a f,. Such an 
expression may be found in [fjj: 

a\(X, z) = 5 a b - ^R a cbd (z)X c X d + . . . , (4.4) 
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Here, R a c bd(z) is the Riemann tensor at the origin. Us- 
ing this in (|2.9p . we find that the flattened lapse has an 
expansion of the form 

N ~ 1 + hiR 1 + -(Roioj + Roikjv k )R i R j . (4.5) 
This allows (|4.3[) to be reduced to 

(q + A) lim (R 1 Ir) (hi + n/h.,) = . (4.6) 

r— >0 

There are two ways to satisfy this equation: either 
q = A = (i.e. qij) = 0), or 



d{hi) 
ds 



= hi + Qi J hj = 



(4.7) 



(|4.3[) isn't the only condition that must be satisfied in 
order for B to remain constant. One must also have 



(q + A) lim 4" [(-RVr 2 ) 3/ In TV 
i — >o os L 







(4.8) 



If <7?A = 0, this is trivially satisfied. Otherwise (|4.7p must 
hold, in which case (14.81) reduces to 



5 



070.7 



R 



0(I\K\J)' 1 



= 



(4.9) 



Roughly speaking, (|4.7|) and (|4.9[) imply that the inertial 
and tidal forces felt by a center-of-mass observer must 
remain constant. 

This leads to a very important conclusion: physically 
acceptable charged currents cannot be dynamically rigid 
with respect to arbitrary reference frames. Unless the 
system is particularly simple, no charged rigid current 
can exist. Although the structure of J a will (usually) 
affect n a and R a bcd{z), it seems very unlikely that quasi- 
rigidity alone could conspire to ensure that they satisfy 
(|4.7j) and (|4.9j) except in a few special cases. 

One consequence of this is that the only quasirigid par- 
ticle fluxes that could exist with generic applied forces 
must contain zero total particles. Nontrivial particle 
fluxes therefore cannot be rigid except for special types 
of motion. 

Interesting rigid currents of electric charge always ex- 
ist, however. This is because although q must still vanish 
in general, the charge density may change sign and J a 
needn't be future-directed timelike. Characterizing such 
vector fields is relatively simple. We start by making 
sure that B is independent of s. Again, setting qip = is 
sufficient. Then 



-B/N 



(4.10) 



It is obvious from (|3.4[) that there always exists a large 
class of B's for which this is possible. Note, however, that 
they all vanish at the edge of the body. This implies that 
1Z cannot change (except perhaps if B = (f = 0). It is 
easy to verify that po also remains constant. 

Besides requiring that B = 0, quasirigidity also de- 
mands that d(A I )/ds. Suppose again that qip vanishes, 



so (p satisfies ([3~T0]) . According to ([2726]) . |2~27|) . and 
(|3.ip . the evolution equation for the three-current never 
become singular. There therefore don't appear to be any 
particularly interesting restrictions arising from the con- 
stancy of A 1 . 

Note however that if 1Z = was not already ensured 
by B = 0, then it is implied by the time-independence of 
A 1 . Although rigidity requires that the boundary of the 
object and its central charge density stay constant in the 
pseudo-Fermi coordinate system we have been using, p 
and j 1 will still depend on s in most cases. Their evolu- 
tion equations are rather lengthy, but easily derived from 
this discussion (supplemented with some straightforward 
bitensor identities). 

To reiterate, we have shown that charged quasirigid 
currents often do not exist. As with Born's original con- 
cept of relativistic rigidity, the idea of unchanging mul- 
tipole moments has significant limitations. Of course, 
the definitions we've adopted are rather strict. Although 
(|4.7p and (|4. 9[) can rarely be satisfied exactly, their left- 
hand sides will often be very small in realistic physical 
systems. When this is true, currents may generically ex- 
ist that are in some sense "nearly rigid." But even when 
viewed from afar, the mechanical properties of such ob- 
jects are likely to be quite different from those of "equiv- 
alent" rigid bodies (which must contain singularities). In 
particular, the computation of gravitational or electro- 
magnetic self-forces in a mechanical system would likely 
be problematic. 

Our definition of dynamical rigidity could also be ar- 
gued to be overly strict in that it requires all of the cur- 
rent's multipole moments to be fixed. This has the ad- 
vantage of producing unique evolution equations, but it 
isn't very well-motivated physically. It would usually be 
more realistic to suppose that only the first few moments 
are constant (or approximately constant). This is closer 
in spirit to the type of quasirigidity discussed in (for ex- 
ample) fl9l |. which was shown to have many desirable 
properties. There does not appear to be any existence 
problem for objects with only "low order" quasirigidity. 
But since the evolution of the higher moments is not 
constrained by such a scheme, the motion of the current 
cannot be unique. 

It should also be stressed that the definition of dynam- 
ical rigidity adopted here is different from the one given 
in [8j, IpJ] . There, the constants of motion were multipole 
moments of the stress-energy tensor rather than moments 
of a relevant current vector. This condition appears more 
restrictive, so we conjecture that a similar non-existence 
result can also be proven for that case. It would then 
follow that the total internal energy defined in [8J] could 
never be constant except in particularly simple systems. 



V. COMPARABLE CURRENTS 

It is often useful to discuss the properties of objects 
in a highly dynamical spacetime with language borrowed 
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from something simpler (such as Kerr or Minkowski). For 
example, one might want to say that a neutron star in 
some complicated numerical simulation was essentially 
the same star that had already been studied in isolation. 
Although defining "same" in this sentence may be quite 
tricky, there are several ways to proceed depending on 
the focus of the problem and the amount of effort one 
is willing to expend on it. This section will introduce 
one method which is particularly simple from a mathe- 
matical perspective, and then study some of its physical 
consequences. 

We define two currents (possibly in different space- 
times) to be "identical" if their multipolc moments evalu- 
ated in some pair of frames are identical. Distinguishing 
quantities related to the second current with underlines, 
we also require these frames to satisfy n J (s) = n J (s), 
v 1 = v 1 , and Q 1 j(s) = Q 1 j(s) for all s. The two currents 
should also have identical charges: q — q. These condi- 
tions automatically imply that the frame components of 
the monopole moments of both currents are identical. 

Equality of the higher moments implies that B(R; s) = 
B(R;s) and ^(R;s) = ^(R;s). Finding the conse- 
quences of these relations will use arguments very similar 
to those given in the previous section. 

As a consequence of both £>'s being equal, (|2.15p and 
(f2~2"0)) imply that 

BR'di In N/N = {q-p) A/A - (q - p) A/A . (5.1) 

In physically reasonable systems, the right-hand side of 
this equation is clearly nonsingular. But B generally in- 
volves an o(r -3 ) singularity. This divergence is softened 
by the fact that R 1 S/lnA / 'A decreases at least as fast 
as r 2 near the origin. But this is not enough to allow B 
to be chosen arbitrarily (using only the constraints dis- 
cussed thus far). In order to remain consistent, it must 
be true that 



(q + A) lim (R 1 /r 2 ) di In N/N 



= 



Using (|4.5[) . this is equivalent to either qip = or 



Roioj + Rq(i\k\j)V 



K 



Rr, 



.0/0,7 + Ro(I\K\J) v 



K 



(5.2) 



(5.3) 



As before, these Riemann tensors are evaluated at the 
respective origins of each frame. This equation essentially 
states that physically reasonable charged currents can 
only be comparable if observers attached to the center of 
each frame would experience identical tidal forces. If two 
charged currents were "comparable" (as defined above), 
but did not satisfy (|5.3p , at least one must be singular. 

We may now construct comparable currents satisfying 
qtp = 0. In this case, 



(jV/A>) <p . 



(5.4) 



It follows that po = P and 1Z = TZ (at least if B ^ 0). 

Continuing, A 1 must be equal to A 1 if two current are 
to be comparable. Suppose again that q4> = 0. If B 



also vanishes, the two current potentials would become 
equal: J 1 = J_ J ■ This would imply that the frame com- 
ponents of the two three-currents would be be directly 
proportional to one another: 



J 



AN~ 
~AN~ 



J 



(5.5) 



In any case, the situation here is similar to the one 
obtained in the previous section. The idea of having two 
objects where all of the relevant moments are identical 
can only be generic if they both carry zero charge. 



VI. DISCUSSION 

We have shown how to construct a vector field J a 
satisfying (|1.2p from a set of multipole moments defined 
with respect to a given reference frame. In so doing, we 
have implicitly found restrictions on the moments which 
ensure that the resulting current will be smooth and 
bounded. This effectively completes the program Dixon 
initiated in 0, [l(| to study the multipole moments of 
nonsingular conserved vector fields in curved spacetime. 

An immediate application of these constructions is a 
method to easily construct currents with a given total 
charge q. This is described in detail in Sec. IIII1 where a 
few free functions with relatively intuitive interpretations 
are used to build any smooth J a . 

In combination with the equations of motion derived 
in one could (for example) derive the equations of 
motion for all charged particles with a given charge in 
some approximation scheme. This was carried out in Q 
in flat spacetime in part to determine precisely when the 
Lorentz-Dirac equation is valid. The results derived in 
this paper would then allow similar methods to be used 
to study the validity of its generalization, the Dewitt- 
Brehme equation 

The methods used here may also be seen as a pro- 
totype for a similar analysis of the multipole moments 
of the stress-energy tensor (as defined in [J, ||). This 
would allow one to directly construct all bodies having 
prescribed linear and angular momenta at a given center- 
of-mass position. 
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APPENDIX A: BITENSORS 

This appendix briefly reviews the theory of bitensors. 
The detailed properties of such objects have been dis- 
cussed in detail by several authors [1, H, [2(| Hl| , so we 
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shall include only a few relevant definitions and identi- 
ties. 

Our discussion will focus on the world function bis- 
calar a. Given two points x and y in a spacetime, it is 
first assumed that there exists exactly one geodesic con- 
necting them. In this paper, these points will always be 
spacelike-separated, so this is a rather mild restriction. 
Regardless, a(x, y) is defined to be one half of the square 
of the geodesic distance between its arguments. 

The derivatives of a arise in many expressions, so it 
is standard to define a special notation for them. In 
particular, let 



a a (x, z) 
a a '(x, z) 



V Q cr(x,z) , 
V Q 'cr(a;, z) 



(Al) 
(A2) 



Here, we have assumed that z 6 Z 1 and used the 
previously-mentioned convention that indices referring to 
points on Z are not primed. This notation trivially gen- 
eralizes to expressions involving any number of deriva- 
tives. For example, <r aba i c = S7 C a'ba&- Derivatives at x 
always commute with derivatives at z, although the or- 
der of derivatives at a particular point cannot generally 
be interchanged (as is usual for covariant derivatives). 

It may be shown that these vectors are tangent to the 
geodesic connecting x and z, and that their magnitudes 
are just the squared geodesic distance [tl \§L I20L [2ll] : 



= 2a 



(A3) 



This allows one to naturally interpret a a and <j a > as gen- 
eralized displacement vectors between x and z. 

Differentiating (|A3|) an appropriate number of times 
also allows one to generate a number of useful identities. 
The simplest of these are o a a a b = &b, CaC a a' = <?V, etc. 
It follows that a a b = 6 b in flat spacetime, or in general 
when x = z. Similarly, a a a i — —6^, in these cases. 

The matrix inverse of —u a a ' appears a number of times 
in Dixon's formalism, so it is convenient to define 



H a a (x,z) := [-a a al {x,z)\ 



(A4) 



This object also appears naturally in solutions of the 
geodesic deviation (Jacobi) and Killing transport equa- 
tions 0, [1[ . Its derivatives have the form 



V v H a a = H a b H< 



' b 
a.V b'c' 



(A5) 



Another bitensor with similar properties is the parallel 
propagator g a a (x,z). This does exactly what its name 
implies: parallel propagates vectors from z to x along the 
geodesic connecting them. The properties of this object 
may be found in jy, [23]. For our purposes, note that 
g a 'a(x,z) = g aa i(z,x). From its physical interpretation, 
it also follows that 



9 a9 V - o b , 

_a a' x a 

g a'9 b - o b 



(A6) 
(A7) 



An important biscalar built from g a a is known as the 
van Vleck determinant A(x, z). It is defined to be 



A{x,z) :=det (~g a> a o a b 
In [6(, it is shown that 



det (g a a 



50) 



so 



A(x,z) = J 44 det (-aV) 
V 9[x) 



(A8) 



(A9) 



(A10) 



The van Vleck determinant plays a fundamental role in 
the focusing of geodesic congruences [y, [22| . But here, it 
will appear as an important factor in the equations giv- 
ing the current vector in terms of quantities constructed 
from the moments. It arises essentially from a coordinate 
transformation performed in the derivation. 

In particular, a number of relevant functions that we 
deal with are defined on the tangent bundle TM re- 
stricted to Z. Occasionally, a base point z will be fixed, 
in which case we can switch freely between considering 
these objects to be functions on TM (really just the tan- 
gent space T Z M) or on the manifold itself (at least in a 
neighborhood of z). This identification will be made via 
the exponential map, so abusing the notation slightly, we 
make identifications of the form f(X, z) — /(exp 2 X) for 
any function / and vector X € T Z M . 

It is then useful to transform integrals over spacetime 
into integrals over T Z M . The inverse of the exponential 
map is just a a . So define 



X a {x,z) := -a a (x,z) 



(All) 



This acts as a "radial vector" between x and z, and from 
(EHI, the Jacob ian of the transformation x — > X is 



dX 



dx 



00) 



A 



(A12) 



Lastly, we note that directly differentiating (|A8|) yields 



V a \nA = -H b b a b ab 



(A13) 



Other properties of the van Vleck determinant may be 
found in alH. 



APPENDIX B: DIXON'S FORMALISM 

In the present context, Dixon's formalism 0, 0, 0, H, 
H, EH defines a set of multipole moments for spatially- 
bounded vectors satisfying (|1.2p . This set is complete 
in the sense that specifying all of the moments allows 
one to completely reconstruct J. If this vector field is 
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interpreted as the electromagnetic current of a charged 
test body, it is possible to show that these moments also 
have the standard interpretation in terms of expansions 
for the electromagnetic force and torque in successively 
higher derivatives of the applied field [1, 0] • 

In most other formalisms of this type (e.g. [13])) a con ~ 
servation law like (|1.2p imposes differential constraints 
on each moment. There are often algebraic couplings be- 
tween the various moments as well. This leads one to a 
system of equations which are impossible to solve with- 
out assuming a priori that only some finite subset of the 
moments are important to problem at hand. While ad- 
equate for many purposes, this removes any possibility 
of reconstructing the current vector from its moments. 
For example, if only a finite number of the moments 
were nonzero, the associated vector field would involve 
the Dirac distribution and its derivatives. Such current 
densities are clearly unphysical. 

Dixon's moments avoid this type of problem by being 
chosen such that there is exactly one differential equation 
imposed on them as a consequence of (|1.2p . This equation 
is trivial, and only states that the total charge does not 
change. The algebraic constraints on the moments are 
also simple solve in closed form, so it becomes possible 
to discuss the entire set of multipolc moments. 

To see how these simplifications work, first note that 
the moments of some source would usually be defined 
by integrating it against a suitable number of "radial 
vectors" referred to some arbitrarily chosen origin. The 
moments should also be allowed to vary in time, which 
is most naturally allowed by requiring that their defining 
integrals be carried out only over spacelike hypersurfaces 
spanning W. Before defining anything, wc therefore need 
a reference frame which foliates the worldtube with hy- 
persurfaces S(s). Each of these must also be given an 
origin z(s) € S(s). Perhaps imagining that these pre- 
ferred points could represent the position of a physical 
observer, we assume that they form a continuous time- 
like worldline Z C W. 

To be more specific, start by fixing the the reference 
worldline Z. One may then specify the foliation by choos- 
ing a unit future-pointing timelike vector field n a (s) on Z. 
E(s) would be defined as the set of all points x £ W such 
that n a (s)a a (x, z(s)) = (see appendix [A] for the defini- 
tion of a a). This naturally generalizes the flat-spacetime 
concept of a hyperplane which is everywhere orthogonal 
to n a . In certain cases, these hypersurfaces will intersect 
each other (inside W), and therefore fail to be a foliation. 
Such instances would rarely occur in practice, however, 
so we shall ignore them. 

Completely specifying this reference frame also re- 
quires fixing the "time" parameter s (up to an additive 
constant). We do this by normalizing the tangent vector 

»«(*(,)) :=i(,):=^£) (Bl) 

such that v a n a = —1. 

One might be inclined to automatically choose n a — 



v a , but more generality is needed for a number of prob- 
lems in relativistic mechanics. In particular, one might 
want to choose the reference frame for the moments to 
be the center-of-mass frame of a body with worldtube W 
[1, S (3 • In that case, n a would point in the direction 
of the body's bulk momentum vector, while v a would be 
the 4-velocity of its center-of-mass line. These two ob- 
jects are physically distinct, and would only coincide in 
special cases. 

Regardless, specifying z(s) and n a (s) is sufficient to fix 
a frame in which the multipolc moments may be defined. 
All we need now is some notion of a "radial vector" be- 
tween x € S(s) and z{s). A particularly natural choice 
for this is the X a (x, z) defined in (|A11|) . which takes the 
form of a vector at z. Accepting this, a multipole mo- 
ment should then be defined by integrating the current 
against a suitable number of these vectors. 

But that cannot be done directly, as J° (x) £ T Z M 
unless x = z. It will be convenient to parameterize all 
possible replacements for the current by a vector-valued 
distribution J a (X,z). We call this the current skele- 
ton (although our definition of it will differ slightly from 
Dixon's). Its precise relation to J a will be given later, 
but for now, we simply assume that it vanishes outside 
of some finite neighborhood of X a = 0. The 2n-pole 
moment of the current is then defined to be 

Q bl - b " a (z) := [ DXX bl ■■■X b "J a (X,z) , (B2) 



where DX := y/ — g(z)d A X . Note the moments are ten- 
sors on Z. The skeleton is defined on the tangent bundle 
TM, and the moments at z are integrals over T Z M . 

It might appear that the moments defined in this way 
depend on the current at all times, but this is not true. 
The support of the skeleton will be found to lie on the 
surface defined by n a X a = 0, which is simply £(s). So 
one of the integrals in (|B2|) effectively drops out, and we 
fall back on the expected three-dimensional definition of 
a multipole moment at some chosen time. 

Note that (|B2|) trivially implies that for all n > 1, 

Q(bl-b„> = Qbl- b n a _ (gg) 

This is the first of three index symmetries that must be 
satisfied by the moments. The others are not so trivial, 
and form the main content of Dixon's formalism. Wc 
simply state that for all n > 2 0, 0, [l(| , 

Ji bl Q bl - b - l[b - a] = , (B4) 
Q^- b ^ = . (B5) 

The monopole moment also has the special form 

Q a (z) := qv a (z) , (B6) 

where q is the total charge defined by (|l.ip . 

We now need to relate J a to J a . In order to do so, 
it will be convenient to think of the current skeleton as 
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a linear functional on the space of all C°° test functions 
with compact support (as is typical in distribution the- 
ory). Given some <E> a (X, z) related to a <fi a (x) in this 
class, we want to know the behavior of 



DX J a (X, z)$ a {X,z) 



(B7) 



for all possible forms of 4> a . This functional should be 
related in some simple way to 

(J,0> := fd 4 x^~g(x) J a '4> al , (B8) 



which determines J a . Note that the (|B7[) depends on z, 
while (J, <f>) does not. It is therefore simplest to link the 
two functionals together by simply integrating out the 
^-dependence: 



<J,0) 



ds / DX J a (X,z)<i> a (X,z) 



(B9) 



Explicitly relating <fi a (x) and $> a (X, z) would now 
complete the definition of the skeleton in terms of J a . 
It is natural to suppose that x — exp z X, but there must 
also be a propagator which maps vectors at z into vec- 
tors at x. The convenient choice for this happens to be 
the bitensor H a a defined by (|A4[) . This is simply the 
identity if the X a are used as coordinates (with a fixed 
z). In general, though, 



**{X,z) :=H a a (ex Vz X,z)(b a ,(exp z X) 



(BIO) 



Taken together, the preceeding equations completely 
define the moments (or equivalently the skeleton) of the 
current vector. To see that they automatically imply 
(|1.2p . choose a scalar test function ip, and then set <f> a ' = 
V a '(y9. With this definition, 



d 



dX a 



</?(exp z X) 



(Bll) 



For any object depending on both X a and z, operator 
d/dX a actually has a very natural interpretation in the 



theory of vector bundles. It is called a vertical covariant 
derivative in and was given the special symbol V* Q . 

Using this notation, (|B11|) may be rewritten as 



Combining this with (IB9|) . 



(B12) 



(B13) 



The divergence on the right-hand of this equation may 
be evaluated by defining the Fourier transform of the 
current skeleton. Let 



(V-J,<p) = ds DXN* a J a ) p 



T[3] a {k,z) := / DXJ a (X,z) 



(B14) 



-ik b X° 



where fc° G T Z M. It then follows immediately from (|B2|) 
that 

T[3] a (k, z) = J2 y -^-k bl ■ ■ ■ fc fc „Q fcl -"" a (z) . (B15) 

n=0 

So the Fourier transform of the current skeleton is a gen- 
erating function for the moments. 

It will be convenient to define a reduced current skele- 
ton J a (X, z) that drops the first (monopole) term in this 
series: T[3] a = J 7 [J] a - Q a . Noting that T[V w A a ] = 
ik a ^F[A] a for any A a (X, z), (|B5j) immediately implies 
that 



V* a J a = 



(B16) 



The divergence of the full current skeleton therefore de- 
pends only on the monopole moment. Using the identity 
T[8{X)\ — 1, we see that 



V* a J a = Q a V*J(X) . 



(B17) 



Substituting this result into 1|B13[) . and using (|B6p finally 
shows that (v a 'J a ,ip^ = for any ip. Hence, (jl.2j) is 
always satisfied, as claimed. 



[1] R. Beig and B. G. Schmidt, Class. Quant. Grav. 20, 889 [10] 

(2003) [11] 

[2] A. I. Harte, Phys. Rev. D 73, 065006 (2006) [12] 

[3] W. G. Dixon, J. Math. Phys. 8, 1591 (1967) [13] 

[4] W. G. Dixon, Phil. Trans. R. Soc. London A277, 59 [14] 

(1974) [15] 

[5] J. M. Sanchez and E. Poisson, | gr-qc/0512lTl~1 (2005) [16] 

[6] E. Poisson, Living Rev. Rel. 7, 6 (2004) [17] 
[7] W. G. Dixon, Proc. R. Soc. London A314, 499 (1970) 

[8] W. G. Dixon in Isolated Gravitating Systems in General [18] 

Relativity, ed. J. Ehlers (1979) 

[9] J. Ehlers and E. Rudolph, Gen. Rel. Grav. 8, 197 (1977) [19] 



W. G. Dixon, Proc. R. Soc. London A319, 509 (1970) 

H. E. Moses, J. Math. Phys. 17, 1821 (1976) 

A. Nisbet, Proc. R. Soc. London A231, 250 (1955) 

E. G. Rowe, J. Phys. A 12, 145 (1979) 

M. Born, Ann. Phys. (Leipzig) 30, 1 (1909) 

N. Rosen, Phys. Rev. 74, 54 (1947) 

G. Salzman and A. H. Taub, Phys. Rev. 95, 1659 (1954) 

H. D. Wahlquist and F. B. Estabrook, J. Math. Phys. 7, 
894 (1966) 

R. Schattner, Gen. Relativ. Gravit. 10, 377 (1979); 10, 
395 (1979) 

M. Barreda and J. Olivert, Int. J. Theor. Phys. 36, 771 



13 



(1997) time (Cambridge University Press, Cambridge 1975) 

[20] J. L. Synge, Relativity: the general theory (North- [22] M. Visser, Phys. Rev. D 47, 2395 (1993) 

Holland, Amsterdam 1966) [23] M. Mathisson, Proc. Cambridge Phil. Soc. 36, 331 (1940) 
[21] F. G. Friedlander, The wave equation on a curved space- 



